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Abstract 

Using analytical methods, we derive and extend previously obtained numerical re- 
sults on the low temperature properties of holographic duals to four-dimensional gauge 
theories at finite density in a nonzero magnetic field. We find a new asymptotically 
AdSs solution representing the system at zero temperature. This solution has vanishing 
entropy density, and the charge density in the bulk is carried entirely by fluxes. The 
dimensionless magnetic field to charge density ratio for these solutions is bounded from 
below, with a quantum critical point appearing at the lower bound. Using matched 
asymptotic expansions, we extract the low temperature thermodynamics of the system. 
Above the critical magnetic field, the low temperature entropy density takes a simple 
form, linear in the temperature, and with a specific heat coefficient diverging at the 
critical point. At the critical magnetic field, we derive the scaling law s ~ T 1 / 3 inferred 
previously from numerical analysis. We also compute the full scaling function describ- 
ing the region near the critical point, and identify the dynamical critical exponent: 
z = 3. 

These solutions are expected to holographically represent boundary theories in 
which strongly interacting fermions are filling up a Fermi sea. They are fully top- 
down constructions in which both the bulk and boundary theories have well known 
embeddings in string theory. 



1 This work was supported in part by NSF grant PHY-07-57702. 



1 Introduction and summary of results 



Understanding the properties of interacting fermions at finite density is a central problem 
in diverse areas of physics, ranging from condensed matter to nuclear and astrophysics. The 
dual gravity description arising in the AdS/CFT correspondence may provide a new window 
on the strongly coupled versions of such systems, and the last few years have seen significant 
effort in this direction. 

The model that has attracted the most attention consists of studying bulk fermions prop- 
agating on a Reissner-Nordstrom black brane background [H 121 El E] • The advantage of this 
model is that it is very simple, yet the fermion correlators display interesting non-Fermi 
liquid behavior. But there are some serious limitations. In general terms, this construc- 
tion does not harness the full power of gauge-gravity duality, as the fermions are described 
"explicitly" in the bulk rather than holographically. In describing boundary fermions we 
would like gravity to do the full path integral for us (at large N), but in this setup the task 
of studying quantum fluctuations of the fermions is left to the user. Another problematic 
feature is that the non-Fermi liquid behavior of the correlators hinges on the existence of a 
near horizon AdS2 factor with finite ground state entropy density. This makes the nature of 
the dual field theory rather obscure, and seemingly quite different from the real world sys- 
tems with which one is hoping to make contact. Other approaches to studying holographic 
fermions, with various pros and cons, include (3, El d E] 

In fact, a promising holographic setup for studying interacting fermions at finite density 
has been sitting right under our noses, in the form of the duality between 5-dimensional 
Einstein-Maxwell-Chern-Simons theory and 4-dimensional supersymmetric gauge theories. 
Simply turning on a finite charge density in this system is inadequate, as it leads to the 
same finite entropy density ground state Reissner-Nordstrom geometry noted above. This is 
presumably somehow a reflection of the fact that the gauge theory contains massless charged 
bosons, which will want to condense at nonzero chemical potential, precluding the formation 
of a Fermi surface. However, turning on a magnetic field removes these problems: in the 
gauge theory the bosonic modes are lifted up in energy above the lighest fermion modes, 
and in the bulk a smooth zero entropy ground state geometry emerges. The presence of the 
magnetic field also provides access to a tunable parameter and, just as in the real world, it 
is interesting to see how the physical properties change as a function of the field strength. 

In previous papers [91 [TUl [11] we have initiated a detailed study of this system, and a 
rich structure has emerged. Perhaps most interestingly, precision numerical analysis revealed 
the existence of a quantum critical point [llj. The results are summarized in Fig. [I] (taken 
from [llj). The system is studied as a function of the dimensionless magnetic field to charge 
density ratio: B = B/p 2 ^. The thermodynamics properties are expressed in terms of the 
dimensionless temperature T and entropy density s. For sufficiently large magnetic field, we 
observed in [H] a linear dependence of the low temperature entropy density, s = c(B)T. As 
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Figure V. Schematic phase diagram illustrating the various behaviors of the entropy density 
versus temperature and magnetic field. The region inside the dashed line is controlled by 
the quantum critical point at (T = 0, B = B c ), and the entropy density can be expressed in 
terms of a sing le scaling function / of (B - B c )/T 2 / 3 . We move around inside this region by 
changing the temperature T and the relevant coupling B — B c . The boundary of the region 
is defined to be where irrelevant operators become important. The yellow region denotes a 
regime where temperature is the largest energy scale, corresponding to the argument of the 
scaling function / being small. Outside the yellow region the low temperature behavior of 
the entropy density, for fixed B, is either constant or linear in T, depending on whether the 
quantum critical point is approached from below or from above B c as T — > 0. 



the magnetic field was decreased, or equivalently the charge density increased, the coefficient 
c(B) was found to diverge at a critical value B = B c . At this critical point a new scaling law 
for the entropy emerged, s ~ T 1 / 3 . Departures from the critical point with respect to both 
temperature and magnetic field could be expressed in terms of a single scaling function, as 
is familiar from the study of both classical and quantum critical phenomena. The numerical 
results suggested that we identify the critical point as having dynamical critical exponent 
z = 3, and possessing a relevant operator of scaling dimension 2. 

Several considerations render this behavior especially appealing. In terms of the connec- 
tion to physical materials, we note that tuning a magnetic field is a common way of locating 
quantum critical points in heavy fermion compounds |T2j , and the appearance of the critical 
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point closely parallels what we observe]^] Namely, away from the critical point one has a 
linear entropy density versus temperature law due to the presence of a Landau-Fermi liquid, 
but the coefficient diverges at the critical point, giving way to a non-linear relation. Another 
helpful aspect is that our behavior is arising in a fully "top-down" setup: the gauge theories 
are explicitly known, and include as one example M = 4 super- Yang-Mills theory. 

On the other hand, since the analysis in [TT] was heavily numerical the basic mecha- 
nism underlying the observed behavior was not entirely clear. Also, numerical results were 
obtained only for a particular value of the bulk Chern-Simons coupling, namely that dic- 
tated by supersymmetry. To both understand and extend our previous results an analytical 
treatment is clearly desirable. This is the basis of the present paper. 



1.1 Results 

We have found a new asymptotically AdSs solution, which is at zero temperature and cor- 
responds to a nonzero magnetic field and charge density in the boundary theory. In a 
convenient set of coordinates the solution takes the formE] 

dv 2 

ds 2 = — + Mdt 2 + 2Ldtdx 3 + e 2V (dxl + dxj) (1.1) 
L 2 

F = bdxi A dx2 + Edr A dt 

Here the magnetic field b is a constant while the functions (L, M, V, E) depend on r only. 
The explicit solution can be written almost entirely analytically: the function V must be 
determined numerically, but it obeys an equation with no dependence on free parameters 
and so is determined once and for all, and then the remaining functions may be solved by 
quadrature in terms of V. In the near horizon limit the functions take the simple form 

2 2k l. 

L = 2br, M = -ar- fc( 9 fc r _ 1} , e 2V = ^ , E = qr k ~ l (1.2) 

where k is the coefficient of the Chern-Simons term in the bulk action, and q is a constant. 
The (r, t,x^) part of this near horizon metric can be recognized as the three-dimensional 
part of the Schrodinger spacetime used in non-relativistic holography [T4"l [T5] . It is also a 
solution of three-dimensional gravity with a gravitational Chern-Simons term, and has been 
studied in this context recently in [16J. 

Viewing this solution as the zero temperature limit of a finite temperature solution, 
we demand that M be negative near the horizon. For k > 1/2, this implies that the 
parameter a must be non-negative. As we show, this inequality translates into a bound on 



2 As discussed in [TT] there is also a close connection with the results described in [T5] . 
By a coordinate transformation the asymptotic metric may be put in canonical AdSs form. 
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the dimensionless magnetic field to charge density ratio measured at the asymptotic AdSs 
boundary. We find B > B c where B c is a fc-dependent number that can be computed in terms 
of integrals involving the V function. This explains the appearance of a critical magnetic 
field in our previous num erical work [TT], and a quantum critical point at B — B c . From 




the form of M(r) in (1.2) it is evident that k = 1/2 is special, and indeed we find that B c 
diverges as k — > 1/2. The quantum critical point is absent for k < 1/2. 

We extract the low temperature thermodynamics by perturbing around the zero temper- 
ature solutions. The appropriate technique employs a matched asymptotic expansion, where 
we match two different perturbations, growing away from the horizon and from the AdSs 
region respectively. This analysis provides formulas for the low temperature entropy density, 
and gives us the dynamical critical exponent since it involves scaling the near horizon time 
and space coordinates as the temperature is scaled. 

Above the critical magnetic field this analysis yields the simple result 

f , B > B c (1.3) 

This result displays all the features that arose in our previous numerical analysis. The linear 
behavior in temperature is characteristic of a Landau-Fermi liquid, here presumably arising 
from the fact that in the field theory we are filling up states in the lowest fermionic Landau 
level. Recall that the specific heat C obeys the same behavior, since C = Tds/dT. In the 
large B limit the result s = 7rT/6 agrees with what we found for the purely magnetic solutions 
in In this limit the near horizon geometry is BTZ xR 2 , and the result corresponds to 
using the Brown- Henneaux central charge [TT] and the Cardy formula. Away from this limit, 
this approach is not applicable, as the near horizon solutions are no longer BTZ black holes 
on account of the nonzero charge density. As we approach B c from above, we see that 
the coefficient of the linear term diverges, as was seen numerically in [TT]. This is directly 
analogous to the observation of a divergence in the specific heat coefficient in real materials 
approaching a field tuned quantum critical point. A finite temperature version of a magnetic 
field induced phase transition was studied holographically in [18], and for zero temperature 
versions involving probe branes see [T9j 



Tuning to the quantum critical point corresponds to setting a = in (1.2). The details 
of the matched asymptotic expansion now change, and the corresponding result for the low 
temperature entropy density is 

s = at 1/3 , a 3 = — B = B C (1.4) 

576fc£3 

This is in accord with the numerical results obtained in [TT]. The scaling analysis shows that 
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the dynamical critical exponent is now z = 30 

We can perturb around the critical point by including both a finite temperature and a 
deviation in the magnetic field away from the critical value. This leads to the scaling form 
for the entropy density 

^ I/3 /5=#) ^ 

where the function / can be found in terms of the solution to a cubic equation, 

Given the result z = 3, this shows that B — B c represents a relevant coupling with scale 
dimension 2 at the fixed point. Again, we have agreement with the numerics in [11]. The 
behavior of other quantities, such as the magnetization, near the critical point may be 



obtained from (1.5) by applying standard thermodynamic relations. 

In other studies of quantum critical points in the AdS/CFT correspondence, nontrivial 
dynamical critical exponents are associated with near horizon Lifschitz metrics [21]. Our 
mechanism is different, and is tightly linked with the fact that our solutions are stationary 
but not static. Entropy density and temperature are measured at the horizon, while time 
and space in the field theory are measured at the AdSs boundary. If the natural time and 
space coordinates at the horizon and at the boundary are nontrivially related, as is the case 
here, new dynamical scaling laws can ensue. This same mechanism can be expected to play 
a role in other AdS/CMT applications. 

A noteworthy feature of our zero temperature solutions is that the charge density mea- 
sured at infinity is carried entirely by fluxes in the bulk; some other recent solutions with 
this property include [221 [21] ■ This is possible by virtue of the Chern-Simons coupling, and 
is in contrast to the Reissner-Nordstrom solution, where all the charge is hidden behind the 
horizon. The fermions in the boundary theory are thus described fully holographically as 
classical bosonic fields in the bulk, and the charge is also accessible since it is outside the 
horizon. It will be interesting to study the dynamical response of these holographic fermions 
by studying perturbations at nonzero frequency and wavelength. 

The remainder of this paper is organized as follows. In section 2, we review the Ansatz, 
reduced field equations, and first integrals, in the presence of a uniform magnetic field and 
electric charge density, exhibit the allowed near-horizon solutions, and review the physical 
parameters of the problem. In section 3, we derive our new charged asymptotically AdS*, 



4 For reasons that will be explained in section 4.3.6 below, this scaling relation holds for all k > 3/4. For 
1/2 < < 3/4 a different scaling relation holds, to be examined elsewhere |23j . 
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solutions, and derive a formula for the critical magnetic field B c as a function of k. In section 
4, we develop a procedure for matching the perturbation theories around the near-horizon 
and asymptotically AdS$ regions, and use this to derive the low temperature formula for 
the specific heat coefficient as a function of B, and the scaling function / near the quantum 
critical point. In section 5, we present a detailed technical discussion of this matching 
procedure, and solve the associated perturbation theories. In Appendix A, we derive the 
most general near-horizon solutions. In Appendix B, we present a detailed discussion of the 
solution of the more involved equations for the asymptotically AdS§ perturbation theory. 

2 Field equations and simple analytic solutions 

The action for 5-dimensional Einstein-Maxwell theory, with negative cosmological constant 
and Chern-Simons term, is given 



Sem = -tt^t ! d 5 x^ (r + F MN F MN -^)+ S cs + S hndy (2.1) 



iG-<r/r, / v -'V u L-'-J 

where the Chern-Simons term is 



For the value k = 2/ a/3, the action coincides with the bosonic part of D = 5 minimal gauged 
supergravity. In this paper, however, k will often be kept general, thus allowing for values 
different from the supersymmetric case. We assume k > without loss of generality, since 
a sign reversal of k is equivalent to a parity transformation. Boundary terms in the action 
are required for finiteness of the action and the existence of a well posed variational problem 
iBJ; their explicit form may be found in [10J. We henceforth set L = 1. 

The Bianchi identity is dF = 0, while the field equations are given by, 

= d* F + kF A F 
Rmn = ^9mn + i 3 F PQ F P Qg MN — 2F MP F N P (2.3) 



With k = 2/V3, the action (2.1) is a consistent truncation known to describe all supersym- 



metric compactifications of Type IIB or M-theory to AdSs [271 [28| [29] . This means that 



solutions of (2.3) are guaranteed to be solutions of the full 10 or 11 dimensional field equa- 
tions (although for non-supersymmetric solutions there is no guarantee of stability). It also 
implies that the solutions we find are holographically dual not just to Af = 4 super- Yang- 
Mills, but to the infinite class of supersymmetric field theories dual to these more general 
supersymmetric AdSs compactifications. 



Conventions: R\ UK = 0*1*, - d v T* K + T^T^ - 1^1% - R\ Xu and R = R^ 
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2.1 Ansatz and reduced field equations 

Our general Ansatz for the metric and field strength reads as follows 

dv 2 

ds 2 = — — — + Mdt 2 + 2Ldtdx 3 + Ndxl + e 2V (dx 2 + dx 2 2 ) (2.4) 

L 2 - MN 6 

F = Bdxi A dx 2 + Edr A dt + Pdx 3 A dr 

where L, M, N, V, E, and P are functions of r only, and B is a constant as a consequence of 
the Bianchi identity. To identify event horizons, it is useful to write the metric in the form 

,9 dr 2 L 2 -MiV l9 / L \ 2 w , , 9 , , , 

= L 2 — MN N ^ + N\dx 3 + —dtj +e 2V (dx\ + dx 2 2 ) (2.5) 

This Ansatz differs in its choice of radial coordinate from that used in our previous papers 
P H01III]- It was motivated by previous work on three dimensional gravity [3DJ [3TJ, and its 
advantages will be become clear as we proceed. 

The metric and field strength are form invariant under the following coordinate transfor- 
mations: 

1. SL(2,R): 

M ( a b .) ( M , ad-bc= 1 (2.6) 



. x 3 ) \c d J \x 3 

2. Scale transformations: 



3. r translations: 



4. Rescaling of x\£- 



r ->■ A 2 r , * ->• r* , ^3 -> ^3 (2.7) 



r-^r + r (2.8) 
xi,2 -> ^xi i2 (2.9) 



The reduced field equations are 

Ml ((iV£ + LP)e 2y ) ' + 2kBP = (2.10) 

M2 ((lE + MPje 21 ")'- 2kBE = 

El L" + 2\/'L' + 4(y" + y ,2 )L - 4PP = 

E2 M" + 2V'M' + 4(V r " + V ,2 )M + AE 2 = 

P3 iV" + 2W + 4(V" + V' 2 )N + AP 2 = 

£4 6^"/ + 8(y') 2 / + Wf -g + AB 2 e~ AV = 

CON (V' ff + 7'/' - 6 + - + P V 4y + MP 2 + 2LEP + iVP 2 = 



where we have defined the following SL(2,R) invariant bilinears 

/ = L 2 — MN (2.11) 
g = {L'f-M'N' 

There is some redundancy in this set of equations: by taking the derivative of the constraint 
equation one can show that one linear combination of the equations is satisfied identically. 
Thus, one could omit one of the second order Einstein equations without loss of information. 

2.2 First integrals 

Eliminating E, P and B between the Einstein equations, we find the following simple result, 

(fe 2V )" = 24e 2V (2.12) 
Therefore, we can, in complete generality, solve for / in terms of V as, 

/(r) = 24e~ 2y M f dr' f ' dr"e 2V ^ (2.13) 



Introducing the potentials A and C respectively for E = A' and P = C, we may integrate 
both Maxwell's equations once, 

(NE + LP)e 2V + 2kBC = 

(LE + MP)e 2V — 2kBA = (2.14) 

The integration constants that arise here have been absorbed into the definition of the 
functions A and C. Forming the combinations M x El — L x E2, N x E2 — M x £"3, 



L x E3 — N x El, and using (2.14), we find three first integrals of motion, 



Xe 2V - AkBAC = A 2A = NM' - MN' 

He 2V + AkBA 2 = /i fi = LM' - ML' 

ve 2V + AkBC 2 — Uq v = NL'- LN' (2.15) 

Here, the combinations A,/i, v are the analogues of angular momentum components for 
SL(2,R) and, just as L, M, N do, they transform under the vector representation of SL(2,R). 
They satisfy the following purely kinematic relation, 

iff = ±f9 + 4(A 2 - fjtv) (2.16) 
Finally, Aq, Ho, an d Vq are the constant values of the first integrals. 
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2.3 Asymptotically AdSs solutions 

In these coordinates, the pure AdS 5 solution takes the form 

dr 2 

ds 2 = — + 2r(-dt 2 + dx\ + dx\ + dx\) (2.17) 
F = 

More generally, we will be looking for asymptotically AdSs solutions, for which the leading 
behavior at large r is given by, 

dr 2 

ds 2 ~ -— + CMrdt 2 + CNrdxl + 2cijdxzdt + cyr[dx\ + dx\) (2-18) 
F ~ bdx\ A dx2 H — %dr A dt H — A 

with c 2 L — cmCn = 4, as implied by the field equations. Here, we denote the magnetic field 
b, reserving use of B for the case where Xi^ are scaled to put the metric in canonical form. 

By using the coordinate transformations discussed above, we have some freedom to change 
the asymptotic parameters. For instance, we can use the SL(2, R) and scale transformations 
to set cl,m,n to any desired values obeying c 2 L — CmCn = 4, and also choose an arbitrary 
value for eg, assuming that it is nonzero to begin with. The freedom to translate r can then 
be used to set the location of the event horizon, if present, at r = 0. 

With b = 0, a familiar asymptotically AdSs solution is the electrically charged Reissner- 
Nordstrom solution, with a near horizon AdS2 x R 3 geometry. In this work we are interested 
in solutions with 67^0, corresponding to a gauge theory with a nonzero background magnetic 
field. 



2.4 Near horizon solutions 

With b 7^ 0, it is possible to find exact solutions where V is constant, so that the metric 
takes the form M3 x R 2 . Such solutions can, and will, serve as near horizon geometries for 
full asymptotically AdSs solutions. In fact, it is not hard to find the general solution with 
constant V. As shown in Appendix |A| for k 7^ 1 (recall that we have restricted to k > 0), 
the general solution with constant V is given by, up to coordinate transformationsP 



F = bdx 1 Adx 2 + qr r>k ~ 1 dr Adt (2.19) 

dv^ ( (p 1 \ b 

dg2 = -JT^ ~ I 3 + ar + , / , 1 \ r2vk dt 2 + 4rjbrdtdx 3 + -^(dxf + dxj) 
w z r z \ k{k - 7,77) J V3 



6 The 77 = 1 solution obviously breaks down for k — 1/2; this value of k is one, among several, "special" 
values of the Chern- Simons coupling k. 
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with rj = ±1. This solution appears to be new, although some aspects are familiar. If we set 
q = the solution is identified as AdS3 x R 2 which is the near horizon part of the magnetic 
brane solutions found in The AdS3 part can equally well be taken to be a nonextremal, 
rotating BTZ solution, since these can be obtained from AdS 3 by applying the coordinate 
transformations discussed in section 12.11 Note also that we can use the residual coordinate 
freedom to set b, q and a to convenient values, such as b = \/3, q — 1 (assuming nonzero q), 
and a = 0. 

A simple but important fact is that the a = /3 = q = solution is invariant under 
the following coordinate transformations, which are a combination of scale and SL(2, R) 
transformations , 

11 , 
r -> Ar , t^—t, x 3 -> — x 3 (2.20) 

for any parameter w. If we write t ~ (x 3 ) z we have 



w 

1 — w 



(2.21) 



Eventually, z will be identified with a dynamical critical exponent. Of course, nothing we 
have said so far determines any particular value of w and hence of z. Preferred values 
emerge once we match the near horizon geometry on to a low temperature asymptotically 
AdSs solution. Doing so, we will find that w — 1/2 and w = 3/4 are realized, leading to 
z = 1 and z = 3 for the values of the dynamical critical exponent in the dual field theory. 

For nonzero q, the (r, t, X3) part of the metric coincides with the corresponding three 
dimensional part of the "Schrodinger" spacetimes, proposed in [T5, 15J in connection with 
non- relativistic AdS/CFT0The nonzero q solutions, but still with a = f3 = 0, are invariant 
under 

r^Xr, t^—t, X3 ^^- k x 3 (2.22) 

Perhaps surprisingly, this scale transformation will end up playing no role in our discussion 
of scaling near quantum critical points; rather, as noted above, it is the scale transformations 
(2.20) with w = 1/2 and 3/4 that turn out to be important]^] Indeed, note that since the 



supersymmetric value of k is greater than 1, the dynamical critical exponent corresponding 



to (2.22) would be negative for the supersymmetric case. 



7 In those solutions there are additional spatial directions appearing in the metric as rdyidyi, and associated 
Galilean boost symmetries. 



In fact, for k < 3/4 the scale transformation (2.22) does determine the dynamical critical exponent 



The importance of k = 3/4 can be seen from the discussion in section 4.3.6. 
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Solutions with the same three dimensional metric, but with the time and space coordi- 
nates swapped, have been studied recently as solutions of topologically massive gravity [16] . 
By virtue of their asymptotic symmetry algebra, they were found to possess a chiral Virasoro 
algebra. 

Precisely at k = 1 a new family of solutions emerges. These are the warped AdS3 x R 2 
solutions found in [10], whose three dimensional part was studied previously in the context 
of topologically massive gravity [22] • These solutions will not play any role in what follows, 
but we present them here for completeness. The solutions are 

dv 2 

ds 2 = ; r -2r(q 2 r + 6r )dt 2 + 4brdtdx 3 + dx 2 : + dx 2 + dx 2 . (2.23) 

12r(r + r ) K J s 1 2 v ' 

F = bdxi A dx 2 + qdr A dt 
with q 2 + 2b 2 = 6. 



2.5 Physical parameters 

Given an asymptotically AdS5 black hole solution we will be interested in computing the 
magnetic field, charge density, temperature and entropy density. These quantities can be 
read off from the asymptotic data near the horizon and at infinity. 



For the near horizon behavior it is clearest to work with the metric in the form (2.5). The 
horizon occurs when L 2 — MN = 0, which we assume is satisfied at r = 0. Finite entropy 
density solutions will have N(0) ^ 0. The value of L at the horizon is interpreted as a 
chemical potential conjugate to P%, the momentum along x%. This follows from the fact that 
the existence of a smooth Euclidean section requires that periodic shifts of imaginary time 
be accompanied by a shift of x% proportional to £(0). We will restrict attention to solutions 
in which the only thermodynamical potentials are temperature and a chemical potential for 
the charge density, and hence we demand £(0) = 0. For the solutions we find, L will vanish 
linearly in r. It then follows that at a finite entropy horizon we require M to vanish; M'(0) 
will be proportional to the Hawking temperature. 



We assume a metric with asymptotic behavior (2.18). To read off physical quantities it is 
convenient to perform coordinate transformations to put the asymptotic metric in canonical 
form. The SL(2, R) transformation 

t -> t - —x 3 (2.24) 

cm 

removes the dx^dt cross term from the metric, i.e. sets cl = 0, while preserving the condition 
L(0) = 0. Next, we can perform scale and SL(2,R) transformations to set — cm = cn = 
c v = 2. 
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After applying these transformations, which of course also act on the near horizon part 
of the metric, it is straightforward to work out expressions for the physical parameters. We 
then express the results in terms of the asymptotic coefficients appearing in the general 
solution, with arbitrary cm, c^, .... The relevant formulas are 



B = — (2.25) 



P 



2 \c v 

T ^ si/db (- M '<°»v^ ') 

Here s is the entropy density with Newton's constant scaled out: s = G5S/V0I. The charge 
density p is determined from the boundary current according to the standard AdS/CFT 
dictionary; see [10] for details. 

These physical parameters are dimensionful, and can be rescaled by a coordinate transfor- 
mation that preserves the asymptotic AdSs metric. It is thus only dimensionless quantities 
that are meaningful, and these are defined as 

B = (2.26) 




T 



B 3/2 

T 
W 2 



Here we are defining s and T in a slightly more convenient manner than in [TT] , where in the 
denominator we used the combination B 3 + p 2 ; this will clearly have no effect on the scaling 
properties since we'll always have nonzero B. 

3 Zero temperature asymptotically AdSs solutions 

We now turn to the construction of asymptotically AdSs solutions representing zero temper- 
ature, finite density, matter in a magnetic field. As we will see, such solutions can be found 
by quadrature in terms of a single universal V function; V is "universal" in the sense that it 
has no dependence on the magnetic field, charge density, or Chern-Simons coupling k, and 
hence can be determined (numerically) once and for all. 
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3.1 Pure magnetic solutions 

We begin by recapitulating the solutions obtained in [5], with a magnetic field but with 
vanishing electric charge. We can choose coordinates such that these solutions obey 

M = N = E = P = (3.1) 

Furthermore, by scaling x% 2 we can set 

b=V3 (3.2) 

The Ansatz is then 

dv 2 

ds 2 = — + 2Ldtdx 3 + e 2V (dx\ + dx\) (3.3) 



F = V3dxi A dx 



2 

so that t and x 3 play the role of lightcone coordinates. 

The problem is reduced to solving for L and V, which obey 

el L" + 2V'L' + A(V" + V' 2 )L = (3.4) 

e4 3 L 2 V" + l e- iV (L 2 e w )'V - l -L' 2 + 3e~ 4V = 
2 2 v ' 4 

con e- v (L 2 e v )'V + -L' 2 + 3e~ 4V -6 = 
v ' 4 



From (2.13 ) we obtain L in terms of V. We choose the boundary conditions L(0) = L'(0) = 0, 
corresponding to a zero temperature horizon at r = 0. Then, 

L{r) 2 = 24e" 2y M f dr' f ' dr"e 2V ^ (3.5) 



and the problem is reduced to finding V . 

As in [9] we are interested in solutions that interpolate between AdS3 x R 2 at small r 
and AdSs at large r. At small r the equations admit a solution of the form 

V(r) = v ir a + v 2 r 2a + v 3 r 3a + ■■■ (3.6) 

with a obeying the quadratic equation 

3a 2 + 3(7 - 4 = (3.7) 

We choose the root such that V is finite at r = 0, and so take 



\/57 1 

o = \- - - « .758 (3.8) 
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Using our freedom to perform a scale transformation on (r, t, X3) we set v± = 1, and then 
find for the first few terms at small r: 



V(r) 
L(r) 



_ 3 (2a + l)(3a-2) r2CT 

3a + 1 

2a - 1 

2 -r CT + • 

a + 1 



(3.9) 



12r 1 



Followed out to large r these initial data match on to an asymptotically AdSs solution. We 
have not succeeded in solving for V analytically, but numerical integration is straightforward. 
The resulting functions are shown in Fig. [2j The large r asymptotics are 



V(r) 




L(r) 100 




20 40 60 SO 100' 



20 40 60 !0 100 



Figure 2: Numerical solution of (3.4) for V(r) and L(r) subject to boundary conditions (3.9) 



e 2y « ~ Cv r , c v « 2.797 (3.10) 
L(r) ~ 2r 

The physical interpretation of this solution was explored in [9]. In the dual field theory 
it corresponds to massless fermions in a magnetic field. The modes occupying the lowest 
Landau level give rise at low energies to a D = 1 + 1 dimensional CFT, which accounts 
for the existence of the near horizon AdS3 factor. In the case of M = 4 SYM, the free 
field central charge was compared with the Brown-Henneaux central charge, and found to 
differ by a factor of ^3/4. Exact agreement is not expected since these solutions are non- 
supersymmetric . 
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3.2 Charged solutions 

We now generalize to include a nonzero electric charge density. Rather remarkably, we 
find that the full solution can be solved by quadrature in terms of the same V function as 
appeared in the pure magnetic solution above. This fact allows us to deduce analytically 
all the most interesting physical properties of these solutions. These properties include the 
existence of a critical electric charge density at which the system undergoes a quantum phase 
transition, and expressions for the low temperature thermodynamics at and away from the 
critical point. 

We describe these solutions in coordinates such that N = P = 0. Examining the field 
equations, we see that L and V obey precisely the same equations as in the pure magnetic 
case, and therefore we can simply carry over the results from that solution. To complete the 
solution we need to solve for E and M. 

3.2.1 Solving for E 

Recalling that we are scaling xi i2 to set b = v^3, equation M2 now reads 

(Le 2V E)' = Vl2kE (3.11) 

which we integrate as 

E(r) = -^^explVuk f - ^' ,, ) (3.12) 
v ' L(r)e 2y W \ Joo L(r')e 2V ( r ) J v ' 



ce is an integration constant, while cy is the same number that appeared in (3.10). It will 
be convenient to define the following function, 

in terms of which the gauge potential, defined as E(r) = A'(r) with A(0) = 0, is given by, 



A( r ) = ^ e ^(r) (314) 
V3k 



The asymptotics of if) are found to be, 



In r i f \/3eo 
r — > Wv) ~ ~j= H — 1=^~ m 



12 V12A; \c v c E 
1 

2cyr 



r — > oo ip(r) ~ (3.15) 
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where eo in the last term of the first line arises from the regularized integral, 



e 



exp < Vl2k / dr' 

I J oo 



1 



This integral is convergent and produces the following asymptotics of E(r) 



r -)■ 
r — > oo 



£(r) 
£(r) 



~ eo r 



fc-i 



(3.16) 



(3.17) 



The asymptotics of ip guarantee that A, as defined in (3.14), automatically vanishes at the 
horizon, as long as k > 0. Note that the ratio e /cE depends only on the properties of the 
purely magnetic T = solution. 



3.2.2 Solving for M 

The remaining equation E2 for M may be integrated by noticing that it is a linear equation 
in M, whose homogeneous part coincides with the equation for L. Thus, we know that we 
can integrate by quadrature. We set 

M(r) = L{r)ip{r) (3.18) 

and obtain the following equation for <p, 

(LVV)' = ~4E 2 Le 2V (3.19) 

In terms of ip, we have 

(LVV)' = -1644^'e 4 ^ (3.20) 
It is straightforward to integrate this equation, and we find, 

where (3 is an integration constant. Integrating once more, we find, 



^=-vr2 + ^Lwh^-^L dr 'w¥^ (3 - 22) 



with a another integration constant. 
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It is now straightforward to evaluate the asymptotics of M, and we find, 



r -> 



r — > oo 



where 



a — a = lQc v c E J(k) 



M(r) fi-ar 



e 2 r 2k 

k(k - \) 



a 



M(r)~-— r 



J{k) 



dr'- 



3 4 v / 3A;^(r') 



2k Jo L(r') 2 e 2V (r') 
The integral J(k) is convergent for all k > 1/2 and positive. 



(3.23) 



(3.24) 



3.3 Interpretation and emergence of the critical magnetic field 

We are interested in solutions that represent zero temperature limits of black hole solutions, 
with a horizon at r = 0. We therefore demand that M vanish at r = 0, which requires that 



we set (3 = On From the expression in (3.22) it is then apparent that 



a < M(r) 
12 ~ L(r) 



< 



a 



(3.25) 



uniformly throughout < r < oo. Since L(r) is positive, M(r) will be negative definite as 
long as a is non-negative. Finite temperature black hole solutions will have M(0) = with 
M'(0) < 0, the Hawking temperature being proportional to — M'(0). Thus, we need a > 
in order for our solution to represent the zero temperature limit of such black holes. In view 



of (3.24), we thus require 



a > l6c 2 v c 2 E J(k) 



(3.26) 



This translates into a bound on the dimensionless magnetic field B. From (2.25), (2.26) 



and (3.23), we have 



B 



3X 1 / 3 n 
4 



1/3 



2/3 
CvCe 



(3.27) 



9 In fact, the /? ^ solutions also have an interpretation as the zero temperature limit of black hole 
solutions, but they necessarily involve a non-zero chemical potential conjugate to momentum along X3. This 
follows since for nonzero j3 we can always shift r such that M(0) = 0, but then L(0) 7^ 0, which is equivalent 
to having such a chemical potential. It would be interesting to explore these solutions further (the neutral 
version is studied in Appendix D of [10]). but here we exclude them since we are working at vanishing P 3 
chemical potential. 
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and so that the bound (3.26) translates to B > B c with the critical field given as a function 
of the universal number cy, and the Chern-Simons coupling k only, 



Numerical integration yields the values shown in Table 1 for B c for various values of k. As 



k 


B c 


.50001 


35.4050648722 


.6 


1.44700934549 


.75 


.916107730288 


1 


.600520361557 


2/V3 


.499424265324 


2 


.264993652464 


5 


.101365592402 


10 


.050219317885 



Table 1: The critical value of the dimensionless magnetic field to charge density ratio for 
selected values of k. 

can be seen from the integral representation, B c diverges at k = 1/2 and goes to zero at 
large k. For the supersymmetric value, k = 2/a/3, we recover the critical value that arose in 
the numerical studies conducted in The interpretation of this critical field in terms of 
the dual gauge theory is of course an interesting question, on which we will comment in the 
discussion section. 

Before turning to finite temperature, we pause to note a potentially confusing point 
regarding these zero temperature solutions. After setting b = \/3 we presented the solution 
in terms of two integration constants, ce and a, and showed how they determined the value 
of B c . It is easy to see however that, up to coordinate transformations, solutions with 
different values of (cE,a) are in fact equivalent. In particular, the SL(2,R) transformation 
£3 — > £3 + ct can be used to shift a, and a rescaling of t and X3 will rescale eg. Even 
though these solutions are coordinate equivalent, it is appropriate to treat them as physically 
distinct when one regards them as zero temperature limits of finite temperature solutions. 
At finite temperature, and in particular when N(0) 7^ 0, solutions with different B c are not 
related by an allowed coordinate transformation: we are not allowed to perform the SL(2, R) 
transformation £3 — > £3 + ct due to the condition that L vanish at the horizon, which in turn 
is mandated by the absence of a chemical potential for P3. In order to have a continuous 
zero temperature limit, we need to keep the full family of zero temperature solutions. We 



19 



will return to this topic in section 4.3.6 where we discuss the approach to criticality in the 
parameters (5, q). 



4 Low temperature thermodynamics 

The solutions constructed in the last section carry nonzero charge density and magnetic field, 
but have vanishing temperature and entropy density. We now want to heat them up and 
study the low temperature behavior of the entropy density. 



4.1 Matched asymptotic expansions 

Since analytic solutions at arbitrary T are not available, we need to proceed perturbatively 
in small T. However, straightforward perturbation theory around the zero temperature 
solutions will not work, as the perturbations diverge at the horizon; at the horizon, the 
change from zero to finite temperature is not a small perturbation. Instead, we need to 
employ a matched asymptotic expansion. The basic idea is well illustrated by a simple 
example. 

Below is the metric for finite temperature D3-branes in asymptotically flat space: 

ds 2 = H-^-fdt 2 + dx i dx l ) + H^f~ 1 dr 2 + H 1/2 r 2 dn 2 5 (4.1) 

r4 4 



To make the analogy with our problem, suppose we only knew the solution (4.1) in two 
limiting cases: 1) r + = 0, which is the zero temperature solution; 2) rCL, corresponding 
to omitting the 1 in H, which yields the AdS5 Schwarzschild solution times the 5-sphere. 
Using perturbation theory, how could we construct the full solution in the regime r + C I 
(the low temperature regime) and extract the low temperature thermodynamics? Note that 
from the knowledge of the AdS Schwarzschild solution one can of course work out its entropy 
density and temperature, but to import this result to the asymptotically flat solution one 
needs additional information, since one does not know a priori how the AdS Schwarzschild 
time coordinate matches onto the time coordinate in the asymptotically flat region. They 
could differ by a scale transformation, which would rescale the temperature; this is a key 
point for understanding the low temperature behavior of the solutions studied in this paper 



10 



In the matched asymptotic expansion approach we consider two different perturbation 
problems. In the first, we perturb around the AdS Schwarzschild solution, requiring smooth- 
ness at the horizon, which leads to perturbations that grow with r. In the second, we perturb 



10 For the D3-brane example, due to the boost invariance of the zero temperature solution there is actually 
no ambiguity in determining s(T) in this manner, but since boost invariance is absent for our solutions, let 
us ignore this fact to maintain a faithful analogy. 
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around the zero temperature solution, maintaining asymptotic flatness, which gives pertur- 
bations that grow at small r. If r + <C L, there is a parametrically large overlap region, 
r + ^ r ^ L where both expansions are valid. In this region we can match up the free 
parameters appearing in the perturbations, and so connect near horizon data with that at 
infinity. For the D3-brane example, the perturbed solution in the overlap region is given by 

ds 2 = ds 2 AdS5xS5 - ^('dt 2 + dx 2 ) + ^dt 2 + ^- 2 dr 2 + ^dr 2 (4.2) 

illustrating the presence of perturbations that grow at small and large r. 

Our problem is similar, with a few extra complications. First, we have two independent 
perturbations, corresponding to changing T and B. Second, our zero temperature solutions 
are not given by a closed form expression. But it is still possible to use this method to 
determine the low temperature thermodynamics. 



4.2 Summary of matching computations 

We defer a detailed solution of the perturbations problem to the next section. Here we 
just summarize the basic results that are needed to determine the thermodynamics. The 
first part of the problem consists of expanding around a BTZ solution, and the second in 
expanding around the zero temperature asymptotically AdSs solutions constructed above. 
These expansions are matched in an intermediate region where both solutions approach 
AdS3 x R 2 plus small perturbations. This AdS3 x R 2 region is the analog of AdSs x S 5 in 
the D3-brane example. Linearizing around AdS3 x R 2 yields the following result 

E(r) = e r fc_1 
P(r) = Po r ~ fc 1 

L(r) = l + 2br-AVz(^^^v + r° +1 -AVz(^^^v_r-° 

M{r) = m + m^r 
N(r) = no + ri\r 

V(r) = v + r a + v-r-"- 1 (4.3) 

where all free integration constants are indicated and b = \/3. The problem consists of relat- 
ing these integration constants to the parameters that appear in the BTZ and asymptotically 
AdS 5 perturbation problems. 
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4.2.1 Perturbations around BTZ 



An exact solution to the field equations is given by BTZxi? 2 with magnetic flux, which we 
write as 

ds 2 = mrdt 2 + 4y3rdtdx3 + ndxl + dx 2 + dx 2 , (4.4) 

12r 2 + mnr 6 1 1 v ' 

F = V3dxi A dx2 
i.e. b = -\A3 and, 

L = 2V3r , M = -mr , N = n , V = E = P = (4.5) 

In order for this solution to match onto a charged asymptotically AdSs solution, we need to 
add a two parameter family of perturbations. One parameter corresponds to adding electric 
charge, and the second to inducing a flow towards AdS 5 . The perturbation parameters are 
q and t> , given by the values of E and V at the horizonp] 

E(0) = q V(0) = v (4.6) 

A smooth perturbation with these boundary conditions can be found analytically in terms 
of hypergeometric functions. For present purposes we just need results for E and V. At 
large r, their asymptotics match onto (4.3) with parameters 

(mn) K 1 (mn) a 

and where 

i2 fc - 1 r(2fc) i2-r(i + 2 ff ) 

Ek - kv(kf ' K - r(i + a)2 (4 ' 8) 

Here, the coefficient j4„ depends only on the k, but not on q, m or n. Its construction will 



be given in (5.25) of section 5.3 



We now match this result to the small r asymptotics of the asymptotically AdSs solution. 
Here we note that eo and v + are nonzero already in the unperturbed solution, and so to 
leading order we need not concern ourselves with the corrections. For E we already defined 
eo to yield the small r asymptotics (see (3.17)) , and for V we chose v + = 1 (see (3.9)). 
Hence the matching yields 

E ' JI e , 7 \-(v + nq 2 A g ) = l (4.9) 



(mn) k 1 ' (mn)°~ 



11 Note that the normalization for the electric charge at the horizon q used in [TT] was chosen with n = 1, 
and differs by a factor of yfn from the normalization used here. 
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4.2.2 Perturbations around asymptotically AdSs solution 

Now consider perturbations around the zero temperature, asymptotically AdS 5 solution. We 
demand that these perturbations preserve the conformal boundary metric, meaning that they 
leave the asymptotic metric constants cl,m,n,v unchanged. With these parameters fixed, and 



with the magnetic field fixed at b = a/3, it is apparent from (2.25) that any change in B 
occurs via a change vo.Ce- 

Up to coordinate transformations, there is a two parameter family of such perturbations, 
corresponding to changing T and B; we denote the corresponding perturbation parameters 
as €t and e B . After turning on these parameters, we extract the small r asymptotics and 



read off the parameters in (4.3). Here we only need the results for M and N, which are 



m = , m 1 = - (a + e T ) , n = C T e T + e B , n x = (4-10) 

We also need to know the change in the large r asymptotics of the electric field, to keep 
track of how the charge density changes: 

r -)> oo E ~ -| , c E = c E0 + +C B e B (4.11) 

where we now denote the electric field coefficient appearing in the unperturbed solution by 
ceo- The coefficients CV and C B depend on the Chern-Simons coupling k, and their values 
may be extracted from the explicit form of the perturbations around the asymptotically 
AdS^ solution. They will be computed analytically in terms of these data in section 5. 

Since M and iV are already nonzero in the unperturbed BTZ solution we can match to 
their values, neglecting the higher order corrections. This yields 

a + €t = in , Ct^t + £b = n (4-12) 
where we recall that a — a = 16cyC 2 E0 J(k). 

4.3 Thermodynamics and scaling 

We now have all the information we need to read off the low temperature thermodynamics. 
The results are different depending on whether we are at or above the critical magnetic field, 
and so we treat these cases separately. 

4.3.1 Low temperature thermodynamics for B > B c 

We first consider sitting at fixed B > B c and lowering the temperature to zero. To stay at 



fixed 5 we set = 0, and then (4.12) yields 

C T t T = n (4.13) 
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To go to zero temperature we need to take er —> 0. The condition B > B c implies a > 0, so 

(4.14) 



then (4.12) gives that m should asymptote to a finite value 

m = a = a — 16c v c 2 E0 J(k) 



From (3.27) and (3.28) this yields 



a_ B 3 
m ~ B 3 — B 3 



(4.15) 



Now, using cm = —a/\/3 (which follows from (3.23)), b = v^3, and the formulas in (2.25) 



(2.26), it is straightforward to evaluate the following results as €t — > 0: 

~ _ 1 m^fcYh 



s = —y/c v na 



47T 



(4.16) 



Eliminating m and a in favor of T and s automatically cancels the dependence on n, and 
we find the following expression for the leading low temperature behavior of the entropy, 



g= W B 3 \ f 
6 \B 3 -B 3 



(4.17) 



as announced in (1.3) of the introduction. For B near and larger than B c , this formula may 



be approximated by 



18 \B-Bj 



(4.18) 



which nicely reproduces the numerical result C3 = 0.045 obtained in formula (3.9) of [TTj 
(where k = 2/Vo), after taking into account the change in the normalization of s and T 
between this paper and [TT] . The exact correspondence is 



C3 



TlBi 



18(1 + B 3 ) 1 / 3 



(4.19) 



giving approximately C3 = 0.04186. 



4.3.2 Low temperature thermodynamics at B = B c 

Next, we examine the low temperature entropy at fixed B = B c , which means a = 0, or 

4 



a = 1644o J(k) = -c 6 v c E oBt 



(4.20) 
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We again set e B = in order to hold ceo fixed, and thus (4.13) holds. From (4.12) we also 
have €t = m, and therefore 



C T m = n (4.21) 

We then find 

s = at 1/3 (4.22) 

with the numerical coefficient given by 

(4.23) 

Since Ctc eo is a universal number, which depends only on k and on the properties of the 
T = purely magnetic solution, so does the coefficient a. In particular, it is independent of 
temperature, magnetic field and charge density. 

4.3.3 Scaling function 

Turning on both ex and e B allows us to explore a two-dimensional region around the critical 
point, corresponding to changing the temperature and magnetic field. Using the matching 
relations, the change in magnetic field to first order in can be expressed as 

B - B c 2C B 2C B , o /i ^ 

— g — = - o e s = o — (^m - n) (4.24) 

B c 3 ceo 3 c^o 

Combining this with 

~ 1 mJcyn 

T = *v (4.25) 

47r A/a 

lets us write the following cubic equation for n 

„• + ^M^n' + f ^^AV „ - (*0!^f» = (4.26) 

C B B r \2C B B r c v 



Here a takes its value at the critical point, a = 16c v c E0 J(k). In (5.52) of section 5, we will 
obtain the following result for C B 

C B = lQU(k) 2 c 2 v c 3 E0 = ^kc 3 v c E0 J(k)B 3 c (4.27) 
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The entropy density is given by s = y/cvna/24, as in (4.16 ). Combining these facts, together 
with the result (4.23), we find that the entropy density can be expressed in the form 




(4.28) 



with the scaling function / obeying 



c • 



For B > B c , and T small, we have x > and large, so that / may be approximated by, 

f(x) * (4-30) 

As we have a lread y obtained an exact expression for the coefficient of 1/ (B — B c ) for the 
ratio s/T in (4.18), we readily derive an expression for a, and thus for Ct in terms of the 



other parameters, 



a 3 = — , C T = 7 27 A , (4.31) 
576kB 3 c %kc 7 v c% Q Bl 



We have already established earlier that these values are in perfect agreement with our 
numerical results of [H] in this regime. For B = B c , we have x = 0, /(0) = a, and we 



recover the critical scaling law of (4.22). 



4.3.4 Dynamical critical exponent 

The change in power law in the s versus T relation at the critical point is a reflection of the 
change in the dynamical critical exponent. This can made precise as follows. As we lower 
the temperature to zero, the near horizon solution is becoming BTZ xR 2 , with parameters 
m and n. Changing the temperature in this regime corresponds to a scale transformation of 



the BTZ solution, under which m and n transform. Under a scale transformation (2.20) we 
have 

m ->• \ l - 2w m , n \ 2w ~ 2 n (4.32) 
For B > B c , taking the temperature to zero at fixed B corresponds to taking n — > at fixed 



m. Thus we have w = 1/2, and then from (2.21) this yields z = 1. On the other hand, 



at the critical point, B = B c , we saw that we should take both m and n to zero with m/n 
fixed. This requires that we take w = 3/4, and so z = 3. 
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4.3.5 Comparison with the numerical results on scaling of 

The value of s/T 1//3 at B = B c was encoded in the parameter c 2 of [UJ, where k was fixed at 
the supersymmetric value k = 2/y/3. The precise relation, taking into account the required 
conversion of normalizations of s and T between this paper and [IT] , is given by 

aB^ , . 

c 2 = £ 4.33 

(1 + 5 c 3 ) 4 /9 

We find the value c 2 = 0.1265, which compares favorably with the numerical outcome c 2 = 
0.11 of P3]. Finally, as B < B c , we find 



/<*> ~ ^ « ~ ^rf (4-34) 



W2k Bl W2k Bl 

This result may be compared with the coefficient c\ of formula (3.7) in [11]. Taking into 
account the conversion of conventions, we find 

ci = -j= , s (4.35) 

4 v / 2T(5 c ) 1 /2(i + 5 c 3 ) 1 /2 

which gives c\ = 0.221 as compared with the numerical value C\ = 0.172 of [11J. 
4.3.6 Approach to zero temperature in the (b, q) plane 

The data used to parametrize the solutions studied numerically in pU] and [H] consisted of 
the charge density at the horizon q, as well as the magnetic field strength b in the horizon 



frame, (A third parameter, denoted by C'(r+), was associated with the momentum P 3 
of the solution, and could be fixed arbitrarily.) In this parametrization, the electrically 
charged, finite temperature, Reissner-Nordstrom solutions correspond to b = 0, while the 
purely magnetic finite temperature solutions of [9] correspond to q = 0, their extremal T = 
limits being reached respectively at the endpoints (b, q) = (0, \/6) and (6, q) = 0). 

For B > B c (and for k > 3/4, as will become clear below), lowering the temperature 
to zero corresponds in terms of (6, q) to approaching the purely magnetic endpoint (6, q) = 
(V3,0). T 



section 



3.3 



Tris circumstance further illuminates the remark made in the last paragraph of 
that the low T thermodynamics, for all values of B > B c , are governed by 
solutions which are all related to one another by coordinate transformations at T = 0. We 
emphasize that even at this endpoint the full asymptotically AdSs solution carries a nonzero 
charge density; what is tending to zero here is the charge at the horizon. 



12 The normalization of the parameter q used in differs from the normalization of the parameter q used 
here by a factor of ^/n, so that \fnqhere = <7[]j] ■ In this subsection only, q will denote the parameter of [TT] . 
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The departure infinitesimally away from the purely magnetic fixed point at (b, q) = 
(v3, 0) caused by turning on a small q is known analytically from our work here. To see 
this, note that the value Vq of the field V at the horizon is equivalent, via rescaling of x\ t 2, 
to a change in b, which is given to leading order by, 5b = \/3 — b = — 2t> . With the present 
definition of q, the matching equations obtained in ( 4.9[ ) read 



-^k=i= e o, ZZZw( w o + g A q) = 1 ( 4 - 36 ) 



mn p x mn 



The parameter eo is kept fixed, since it is related to the value of B, which is being held fixed 
as we lower the temperature. Rearranging and trading vq for Sb, we have 

q = i2_7ra fc -V-2 , Sb = 2A q q 2 - 2(V a )- 1 {mnY (4.37) 



E 



The parameter A q , which will be obtained in section 5.3, is negative when the ^'-integral is 
convergent, namely for k < 1 + a/2. These relations allow us to examine the curve in the 
(b, g)-plane along which the T = limit is obtained while keeping B fixed. 

• For B > B c , the zero temperature limit is obtained by taking n to zero while holding 
m fixed. It is then appropriate to write 

5b = 2A q q 2 - 2V7 1 m^r (—) ^ ? H (4.38) 

The first term dominates for | < k < ^y^; the second for < k. In either case the q — > 
limit is smooth, since both powers of q are positive. Thus, as T — > 0, the B > B c system 
flows to the purely magnetic critical point for all k > 1/2. Numerical study confirms these 
analytical results. 

• For B = B c , the zero temperature limit is obtained by taking both m and n to zero 
while holding the ratio m/n fixed. We write 

5b = 2A q q 2 - 2V- 1 i^j 5 (^-J q 1 ^ (4.39) 

For | + | < k the second term dominates, and thus determines the flow. For | < k < f + § 
the first term dominates. 

For | < k, both terms admit a smooth q — > limit since, just as in the case B > B c , 
both powers of q are positive. Thus, as T — > 0, the B = B c system again flows to the purely 
magnetic critical point, a result again confirmed by numerical study. 

Finally, for | < k < | the power of q in the second term turns negative, indicating 
that the system no longer flows towards the purely magnetic critical point as T — > 0. In 
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fact, numerical study shows that as T — > 0, a finite limiting value for q emerges. Thus, the 
perturbative expansion around the purely magnetic critical point, which was used throughout 
this paper, can no longer be valid for the B = B c and | < k < | system, and the scaling 

law s ~ Ts is not expected to hold. The precise nature of the flows for k < | is presently 
under investigation, and will be discussed elsewhere [23] . 

The critical curve, defined as the set of limiting values in the (b, q) plane at which T = 0, 
may be computed from the above result, in the neighborhood of the pu rely m agnetic critical 
point, where it may be defined as the flow at B = B, 

by 5b q 2 for 3/4 < k < 3/4 + a/2, and by 5b ~ 

found to be in accord with numerical results. 



We deduce from (4.39) that it is given 
. g <r/(4fc-3) when A; > 3/4 + a/2, a result 



5 Matched asymptotic expansion: detailed analysis 

To explore the low temperature thermodynamics in the regime B > B c , we proceed per- 
turbatively in small T, namely T <C \[B~ C , as outlined in section |4j For r ^> T, the effects 
of finite temperature are small and may be treated perturbatively; in this region, we use a 
perturbative solution around the charged, T = 0, asymptotically AdS 5 solution constructed 
in section [3] For r <C \[B~ C , the finite T asymptotically AdS§ solution reduces to BTZ 
xR 2 ; in this region, we use a perturbative solution around finite temperature BTZ. In the 
overlap region T <^ r <^ \[B~ C the effects of temperature are small and a perturbed AdS^ 
solution may be used. It is in this region that both the near-horizon perturbed BTZ solution 
and the asymptotically AdS§ perturbed solutions are both valid, and where their pertur- 
bative solutions may be matched. In this section, we present detailed derivations of the 
perturbative expansions in each one of these regions, and of their matching. This section is 
somewhat technical, and so we note that the main results have already been summarized in 
the preceding section. 



5.1 Near- horizon region: perturbed BTZ solution 



In the near-horizon region, we expand around the BTZ xR solution of (4.4) and (4.5). 



Inspection of the reduced field equations (2.10) instructs us to treat E and P to first order 



but the other fields to second order. To organize this perturbation theory we introduce a 
small parameter e, so that 



E(r) 
P(r) 
L(r) 
M(r) 



eE 1 (r) 
ePi(r) 

2br + e 2 L 1 {r) 
—mr + e 2 Mi{r) 
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N(r) = n + e 2 N l (r) 

V(r) = eVi(r) (5.1) 

where b = v^3, and m and n are the constant parameters of the BTZ solution. The corre- 
sponding field equations are, 



Ml 


= 


nE[ + 2brP[ + 2b(k + 1)P 1 


M2 


= 


2brE[ + 26(1 - k)E 1 - mrP[ - mP x 


El 


= 


(Li + 8brV; - AbVj - AE X P X 


E2 


= 


(M[ - AmrV{ + 2mV l )' + 4^ 


E3 


= 


(N[ + AnVi)' + 4P X 2 


EA' 


= 


6((12r 2 + mnr)V^' - 96Vi + AnE\ + \§brE x P x - AmrP? 


CON 


= 


in 

bL[ + — N[ + (24r + mn)Vj -12]^ + nE\ + AbrE x P x - mrP 2 



We have replaced E4 by E4'= E4 + 4 CON. The boundary conditions at the horizon are, 

Li(0) = Mi(0) = M[(0) = iVi(O) = (5.3) 

The first two may be chosen by SL(2,R) and the requirement that the horizon remain at 
r = 0; the latter two follow from the fact that any non-zero values may be absorbed into the 
parameters m,n of the zero-th order solution. The values of Li(0) and Vi(0) will be turned 
on by the perturbation, and -Pi(O) will be determined by those, as will be shown below. 



5.1.1 Solving Maxwell's equations 

Maxwell's equations Ml and M2 involve only E\ and P\ and may thus be integrated inde- 
pendently from the remaining equations. Throughout, it will often be convenient to use the 
rescaled coordinate z defined by, 

(5.4) 

mn 

We shall use the same notation for a function of r and its associated function of z. Eliminating 
P{ between Ml and M2, it follows that Pi is uniquely determined by E 1: 



Pl - 12 



;i - z)E[ + (k - 1)E 1 



(5.5) 



km 

while Ei satisfies the hypergeometric differential equation, 

z(l - z)d 2 z E l + [c - (a + b + l)z]d z E 1 - abE 1 = (5.6) 
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with a + b = 2, ab = 1 — k 2 , and c = 1. Retaining the solution which is regular at r = 0, 
we impose a boundary condition Pi(0) = q, corresponding to finite charge density at the 
horizon. The solution is then given by the hypergeometric function, 

E 1 (r) = qF(l + k,l-k;l;z) (5.7) 

The functions E\ and P\ are both regular as r -> 0, taking values 

Pi(0) = q Pi(0) = -VUq(k - l)/m (5.8) 



5.1.2 Solving Einstein's equations 

Einstein equation E4' may be solved for V\ in terms of E\ and P\. To do so, it will be useful 
to introduce the composite function A, defined by, 

nq 2 A = — (nE 2 + ^brE x P x - mrP 2 ) (5.9) 
18 ^ ' 

By inspection, it is clear that the function A depends only on z and k. In terms of these 
variables, E4' for V\ may be expressed as an inhomogeneous hypergeometric equation, 

z{l - z)d 2 z V x + (1 - 2z)d z V x + \v x = nq 2 A(z) (5.10) 

The homogeneous part is the hypergeometric equation with a + b — 1, ab — —3/4, and c = 1, 
following the notation of (5.6). As a result, a and b are given by a = 1 + a and b = —a, 
where a was defined in (3.8). There is a unique solution which is regular at r = z = 0, which 
we shall abbreviate by, 

F(z) = F(a + l,-a;l;z) (5.11) 

To construct the inhomogeneous solution, we set V\(z) = F(z)v(z), after which d z v is found 
to satisfy a first order equation, 

d z (z(l - z)F 2 d z v) = nq 2 FA (5.12) 

This equation may be readily integrated, and the resulting general solution which is regular 
as z = r = is given by, 

V^z) = v F(z) + nq 2 F(z) £ dz? _ * jT dz" F(z")A(z") (5.13) 

Here, t>o is an arbitrary integration constant. Since E\ and Pi are regular as z — > 0, both 
integrals are convergent, so that V\{z) is a well-defined function of z with Vi(0) = v . 
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Einstein's equations El, E2, E3 may be solved in terms of Ex,Px,Vx, as follows, 

Lx(r) = 4a L r - 8brVx(r) + 12b f dr'Vxir') + 4 f dr' [ dr"ExPx(r") 

Jo Jo Jo 

rr rr rr' 

Mx(r) = 2vqt + 4mrVi(r) — 6m / dr'Vx(r') — 4 / dr' / dr"E\{r") 

Jo Jo Jo 

Nx(r) = Aa N r - An(Vx{r) - v ) - 4 f dr' f dr"P?(r") 

Jo Jo 



(5.14) 



These solutions fulfill all the boundary conditions of (5.3). The two remaining integration 



constants aL,ctN are related by enforcing the constraint equation CON, and we find, 



nq 



-AboiL — maw 



(5.15) 



In summary, given m, n of the unperturbed BTZ solution, the perturbation theory around 
BTZ is governed entirely by two additional parameters q, Vq (the parameter will be 
determined by further gauge fixing and matching). 

5.2 Overlap region: perturbed AdS^ solution 

In the overlap region, the unperturbed geometry is simply AdS^ x R 2 , whose fields are given 



in (3.3) with V = 0, and corresponds to setting m = n = in the BTZ equations (5.1) and 



(5.2 ). All fields are now expanded to first order only, and we find the leading order solutions, 



E(r 
P(r 

L(r 

M(r 
N(r 
V(r 



e r 
p r~ 



k-i 



-k-X 



l + 2br - AV3 

mo + ITLxT 

n Q + nir 



2a - 1 
a + 1 



v + r a+1 - AV3 



2a + 3 



a 



v+r + V-T 



-a-X 



(5.16) 



Note that the constraint equation CON fixes the linear term in L. There are 9 integra- 
tion constants, e ,p ,l ,m ,mx,no,nx,v + ,V-; their number precisely matches the expected 
number from 4 second order equations, 2 first order equations and one constraint. 



5.3 Matching expansions in near- horizon and overlap regions 

To carry out the matching between these two regions, we need to isolate, in the large r/ (mn) 
behavior of the perturbative solution around the BTZ solution, those functional dependences 



that coincide with those identified in the overlap region in (5.16). These large r/{mn) 
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asymptotics arise from the asymptotic behavior for large z of the hypergeometric function, 
which may be obtained using the inversion formula, 

F{a;b;c;z) = - a \ {-z)- a F (a; 1 - c + a; 1 - b + a; z' 1 ) (5.17) 



z)- b F (b;l -c + b-l -b + a^z' 1 



V{b)T{c-a) 
r(c)r(o-6) 

r(a)r(c-6) 

which results in the following dominant asymptotics as —z — > oo , 

, \ T(c)T(b — a) , . „ T(c)T(a — b) , . h , 

F ^ 6; c - z > ~ rwl -„; <-'>-• + fUfW<- 2 » < 5 ' 18 > 

Isolating the r fc_1 term in i?i and the r~ fc_1 term in Pi, we find the following relations, 

a(mvV- k F F - (lg)^r(2fc) 

po = q{mn) l+k E_ k (5.19) 

Isolating the r a and r _1_cr terms in V\ proceeds analogously, paying close attention to some 
extra subtleties. We begin with the asymptotics as r/{mn) — » oo of the corresponding 
hypergeometric function, 



F{z) ~ K(mn)-V ff + F_ 1 _ CT (mn) 1+CT r- 1 - <7 K = 1 > n \ \ 2 ' (5.20) 

1(1 + a) z 

The subtlety in evaluating the asymptotics of V% resides in the fact that contributions arise 
due to the inhomogeneous part of the solution in (5.13). To identify these contributions, 
one needs to investigate the asymptotics of the integrals for large —z. The behavior of A is 
readily deduced from that of E\ (given in ( |5.16[ )), and that of Pi, given by 

mk\mn) K 1 

As a result, the —z — > oo asymptotics of A is given by, 

A(z) ~ Acoimn) 2 - 2 ^ 2 *- 2 A M = — — % (5.22) 



/\2fc-l+c 



The asymptotics of the integral over z" in (5.13) as —z — » oo is proportional to (— z') 
For k > 1/2, as we have been assuming throughout, this exponent is positive. The leading 
asymptotics of the z' integral in (5.13) is then governed by the asymptotics of the integrand, 
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namely (—z') 2k 3 ° . For 2k — 2 < cr, the integral converges as z — > oo, and the leading 
asymptotics of V\ is given by 

V x {r) ~ A v V a {mn)- a r <T + A v V^ r7 (mn) 1+a r- 1 - a (5.23) 

where was defined earlier, and Ay is given by, 

A v = v + nq 2 A q (5.24) 

where A q is a quantity that depends only on k. For 2k — 2 < a (or equivalently A; < l + a/2 ~ 
1.38, so that this range clearly includes the supersymmetric value of k), we have 

^1 " rtl-zDFM' L ^"W) (5-25) 



For larger values of k > 1 + a/2, there is still a contribution of the form (5.24), but it must 
now be obtained after subtracting out the leading r — > oo asymptotics of the integral. For 
example, when a < 2k — 2 < er + 1, we have 



A q = lim 

— z— >oo 



/ dz — r — ; — — - / dz"F(z")A(z") — 

Jo z'(l- z')F(z') 2 Jo v ; v ; 



F{z) 



where 

3A r 



(5.26) 



(5.27) 



2(<ok 2 - 9k + 1) 

Every increment in the range of k by 1/2 will require an extra subtraction term. 

We are now ready to match the asymptotics of V\ between the near-horizon and overlap 
regions, and we find, 

v + = A v V a {mn)- CT 

v _ = A v V^ a (mn) 1+a (5.28) 

The ratio V-/v + = (V-\- a /V a ){mn) 1+2a shows that the term proportional to t>_ corresponds 
to a higher order correction in the asymptotically AdS$ region, and may be neglected there, 
so that we effectively have t>_ ~ 0. 

It remains to match the expansions of the functions L\,Mx,Ni . To leading order, we 
have the matching, 

k = 

mo = mi = —m 

n = n m = (5.29) 

The large r/(mn) asymptotics of the functions L\,M\,Ni may be derived in an analogous 
fashion, and will provide higher order corrections to the leading matching of (5.29). The 
corresponding results will not, however, be needed for the thermodynamic questions that we 
are addressing here, and we shall not carry them out. 



34 



5.4 Asymptotically AdS§ region: perturbed T = charged solution 

In the asymptotically AdS$ region, perturbation theory is carried out around the T = 
charged asymptotically AdS 5 solution, which was derived in section |3j Here, we shall denote 
the fields of this solution with subscripts 0. The functions Lq and Vq obey the equations, 

= L£ + 2^ + 4L (< + a/ ') 2 ) 

= L 2 M) 2 + V^Ll)'-6 + \(L f ) 2 + 3e-^ (5.30) 



for both neutral and charged solutions. For the charged solutions, we have Eq = A' with, 

yc E 
kb 



A ( r ) = CyCE0 c 2kH(r) 



a 2 2 . r 4kbip(r') 

W) = -^W-^lr)/^ ^, (5.31) 



*'< r) s / (5 - 32) 



where the function ip is defined by 



as well as N = P = 0. For given Vb, the first equation of (5.30) is a linear second order 
differential equation. Besides it solution Lq, it has a conjugate linearly independent solution, 
which we shall denote by Lq, and normalize by, 

r dr' 

L ^ sL ^L L (rre^ (5 ' 33) 
This function obeys the following asymptotics, 

1 

r — >■ Lq(t) -> 



r -> oo Lg(r) — (5.34) 



12 
1 

4cyr 



This solution will play a key role in the sequel. 



5.4.1 Perturbation equations 

To carry out first order perturbation theory around this solution, we introduce an expansion 
parameter e, and set, 

P(r) = e P x [r) 
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N(r) 


= eJVi(r) 


E(r) 


= £ (r) + e£i(r) 


L(r) 


= L (r)+eLi(r) 


M(r) 


= M (r)+£Mi(r 


V(r) 


= 7 (r)+e7!(r) 


M 


= /o(r) + e/i(r) 



(5.35) 

Maxwell's equations for the perturbation functions are given by, 
Ml = (e 2Vo (£ iVi + LoPi))' + 2ttFi 

M2 = (e 2Vb (2L P ^i + ^o^i + E Q L X + M Pi))' - 2kbE 1 (5.36) 

while Einstein's equations are given by, 

El = L'[ + 2V Q 'L'i + 2^ + 4L X (V " + (K ') 2 ) + 4L (Vi" + 2y Q V/) - 4E P 1 

E2 = Mi + 2V 'M{ + 2V(M' + 4Afi (v o " + (V^') 2 ) + 4M (V" + 2V£V?) + 8£ £i 

£3 = < + 2V^N[ + 4^ (V " + (\/ ') 2 ) 

/V = /{' + AVJ[ + 4^/0 + 2Vl'f + 2V "/i + 4(y ') 2 /i + 8V^V(f 

EA' = (6V£ + 12(V ') 2 ) /1 + (6V» + 24^%') / Q + 6^'/£ + 6^'/J 

-32fcV 4y °Vi + 8L E P 1 + 4^^ (5.37) 

Here, the last equation corresponds to the combination E4' = EA + ACON, while the next- 



to-last equation derives from expanding the f,V equation of (2.12). 



5.4.2 Boundary conditions 

The T / perturbations on the T = asymptotically AdS§ solutions will be required to 
have specific boundary conditions at r = 00, namely that the space-time retain the same 
asymptotic metric as the T = solution. This will require Li(r)/r, M\(r)/r, Nx(r)/r, 
r 2 E\(r), and V\(r) to tend to as r — > 00. At small r, namely T <r <C B c , we will impose 
matching conditions with the perturbations around the T^O BTZ solution, as given in the 



overlap region by (5.16). 



5.4.3 The translation and dilation modes 



Invariance of the equations (2.10) for the T = unperturbed asymptotically AdS$ solution 



under translations and dilations in r guarantee the existence of two independent perturbative 
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solutions to (5.36) and (5.37). They will be denoted with the superscripts t and d respectively, 
and are given by, 



El = E' Ef = rE' + E /2 

L\ = L'q L\ = rL' — L 

M{ = Mq M? = rM' - M 

VI = V ' V? = rVo' 

fi = fo fi = rfo - 2/o (5.38) 

with P\ = Pf = N\ = Nf = 0. As expected, the translation mode obeys the desired 
boundary conditions at r — > 00, but the dilation mode does not, because Vf(r) and r 2 Ef{r) 
do not tend to there. Thus, the dilation mode must be absent altogether. 

5.4.4 Solving and matching for N\ 

The function N\ satisfies the same differential equation as Lq does. As a result, it must be 
a linear combination of L and its conjugate Lg. The boundary conditions at r — > 00 forces 
the function L to be absent, so that we have, 

iVi(r) =n L c (r) (5.39) 

with asymptotics given by, 



12 

r^oo iVi(r) (5.40) 

4cyr 

Matching the r — > asymptotics with those of the overlap region, we find that 

n = -^L ni = (5.41) 

v 12 



since the function Lq does not contain a dependence on r which is linear for small r. 
5.4.5 Solving and matching for Pi 



The function P 1 = C[ may be obtained by solving Ml in (5.36). It is instructive, however 



to solve for C\ from the first integral for A in (2.15). Linearizing this equation gives 



{N X M' Q - N[M )e 2Vo = 2A + (5.42) 
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This equation is easily solved, and we find the general solution, 



~ A Ll n M A 
Gl " ~ n °^ ~ 8fc&Z^ ~ AkbA {bA6) 

The relevant asymptotics of Pi as r — > is the functional dependence r~ l ~ k , resulting from 
the dependence r~ k of C\. The first term above behaves as r k ~ 1 which has no overlap with 
r -i-h f OI > The last two terms do contribute, however, and we find, 

r — > Ci ~ r 

48e 

^ rioa — 46An ( kb\ . . 

r oo d ~ ^— °- 1 + 5.44 

As a result, we find, 

48e 

Poo n a - 46A 

9 Poo i r 1 

r 2 16c(rC E o 



r ->■ oo Pi ~ ^ Poo = -fe 1g ^ 2 (5.45) 



Comparison with the functional behavior in the overlap region, we find that the prefactors 
Po must be the same in both expressions. The value of po in the overlap region has already 
been determined, and found to be given by po = gP_fc(mn) 1+fc . This value provides a higher 
order correction and vanishes to leading order, po ~ 0. Combining with the result from 
(5.45), we must have to this order, 

A = fi Q a/(Ab) (5.46) 

so that 

Poo = -n k C \ a =-n kc E0 J(k) (5.47) 
lQc 2 v c E0 

The large r asymptotics of Pi will play a key role in determining the constant Cb- 



5.5 Calculation of the coefficient Cb 

With the matching information on N\ and Pi in hand so far, it is already possible to compute 
the coefficient Cb which plays a key role in determining the scaling properties near the 
quantum critical point. Since Cb governs the response of the charge density ce to a change 
in the parameter n at the horizon, we need to determine the change in the field Pi caused 
by turning on the perturbations Ni and Pi. We begin by solving for E\. 
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5.5.1 Solving for E\ 



The function E\ may be obtained by solving equation M2 of (5.36), in terms of the functions 



P\,V\^L\. In fact, it is instructive to obtain the corresponding potential A\ from the first 



integral in the second equation of (2.14). Its perturbative expansion is given by, 



L e 2Vo A\ - 2kbA 1 = -e 2Vo (2L E Vi + + M Pi 



(5.48 



Dividing both sides by A Q L e 2V ° and using (3.14), we readily integrate this equation, and 
the general solution is given by, 



A 1 (r) = Ao(r) 



(5.49) 



where a\ is an integration constant which may be absorbed into the zero-th order charge eo 
or ceo- The matching for small r is carried out between the zero-th order asymptotically 
AdS 5 solution E and E{r) 



e r 



k-1 



in the overlap region (5.16), so that ceo and e are 



related by the zero-th order formula (3.16). The corrections provided by E\ only add higher 



order effects, which may be ignored for our purposes. 

To compute the change in charge density ce, we need to obtain only the leading correction 
to the 1/r asymptotics of A\. Remarkably, to do so, we do not need to know the precise 
solutions for the functions L\ and V\. The only property of Li and V\ that we will need 
here is the fact that they preserve the asymptotic form of the metric, namely as r — > oo, the 
functions rV\ and L\ must have finite limits. 

As a result of the above properties of V\ and L 1; the leading asymptotics of the integrand 
of (5.49) derives entirely from the term in P 1; and is of order 1/r 2 . The integral converges 
at oo, and the asymptotics there is given by, A 1 (r) ~ — ap OQ /(\/l2r), where p^ was defined 



in (5.45). As a result, the asymptotics of E\ is given by 



EAr) 



CEl 



cei 



lQn c v c E0 kJ(ky 



(5.50) 



Thus, the change in the charge density caused by turning on hq is indeed non-trivial. 



5.5.2 Calculating Cb at the critical point 

The coefficient Cb is defined as the change in ce due to turning on iV at the horizon, 

ce = ceo + Cb^b 

N = e B r-+0 (5.51) 

while keeping the temperature fixed at T = 0. Thus, we need to know the asymptotics of 
the field E\ due to turning on n . This is precisely the effect calculated above, and those 
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results may now be applied. The quantity ce — ceo m ay be identified with cei above, and 
e B with n ~ n to leading order, so that we find, 



C 



B 



cei 
n 



kJ(k)acEo = l6kJ(k) 2 CyC E0 



(5.52) 



using the results of (5.50). As was shown in section 4.3.3, once Cb is known, the remaining 



coefficient CV may be calculated indirectly by matching the scaling function near the quan- 
tum critical point onto the low T behavior of the B > B c branch, which is known exactly. In 
the remainder of this section, however, we will solve for the perturbations around the T = 
asymptotically AdS^ solution also for the remaining functions Vi, Li, and Mi in order to 
confirm that the expansion and matching procedure works consistently to leading order also 
for these functions. 



5.6 Solving and matching for Vi, Iq, and Mi 

The calculation of the perturbations V\ and L\ is considerably more complicated than for N\ 
and P\, because their equations do not decouple. In Appendix B, the perturbations of V and 
of the composite f = L 2 — MN are obtained from the perturbations of a common field X, 
via the relations, e 2V = X" and fe 2V = 24X. The zero-th order solution corresponds to the 
L and Vq fields of the T = asymptotically AdS$ solution, and is given by 24X = Lle 2V °, 
since No = 0. The perturbations V\ and L\ are obtained from the perturbations X\ of X, 
and the known perturbation Ni by the relations, 

2V! = X"e~ 2Vo 
2L £i = -2ViLo + 24X 1 e- 2Vo + M A^i (5.53) 

It is shown in Appendix B that the field X\ obeys a third order linear differential equation, 
whose inhomogeneous part is sourced by N\ and P\ . Two linearly independent solutions are 
known to the homogeneous part of this equation, namely the translation and dilation modes 
of X , which we denote by X\ and Xf respectively. They are given by, 

X\ = X' X* = rX' - 2X (5.54) 

The third solution X" to the homogeneous equation, and a particular solution Xf to the full 
inhomogeneous Xi-equation are constructed in Appendix B, so that the general solution for 
Xi, V\ and L\ is given by, 

X 1 = CtXl + CdXf + CnXf + noXf 
V± = CtVl + CdV^ + CnVr + noVl 

L x = CtLl + QLf + CnL^ + noLl (5.55) 
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The translation and dilation modes of V\ and L\ (as well as of the remaining fields) were 
already given in (5.38). The n and p modes are given by 



2V? 
2V 1 P 



{xir 



ny/g-2V 
-2V 



2L L™ = -2V?Ll + 24X 1 n e~ 2Vb 

2LqL\ = -2VfL 2 + 2AXfe~ 2Vo + M N X 



(5.56) 



The large r asymptotics of the X± modes was evaluated in (B.26). The modes t,n, and p 



all respect the asymptotically AdS$ metric, and are thus allowed from the large r boundary 
perspective, while the mode d violates those boundary conditions and must be absent, so 
that Q = 0. As r — > 0, the asymptotic behavior of the remaining modes Xl,X™, Xf was 



evaluated in (B.28), and we derive the following asymptotics for the corresponding V± and 
L\ components, 



,-0—1 



U[ ~ r - a 

L\ ~ r CT+1 



(5.57) 



The V™ behavior as r — > requires it to match on to the r branch of the solution (5.23) 
in the overlap region. As was shown in ( |5.28 ), however, the coefficient t>_ of this branch is 
suppressed by a factor (mn) 1+2ff ~ (mn) 2 ' 5 compared to the v + branch and thus corresponds 
to a higher order effect, which may be neglected. Thus, the solution in the overlap region 
instructs us that the V™ mode must be absent to leading order as well, so that ( n = 0. It 
remains to match the mode V/\ Combining it with the leading contribution from the zero-th 
order asymptotics r CT , we have 



v+ = (v + nq 2 A q )V a (mn)' a = 1 + 2h Q aJ 2 (k) 



(5.56 



where the coefficient J2(k) is defined in (B.29) of Appendix B. Clearly, the second term on 
the rhs is subdominant, and may be dropped in the matching process, yielding the matching 



condition of (4.9). 



5.6.1 Solving and matching for Mi 

The homogeneous part of the only remaining equation E2 is solved by Lq and Lq. The 
general solution is given by, 

Mi(r) = rh L c (r) + rhiL (r) 

-L c (r) £ dr'L e 2Vo (W{M' Q + 4M (K 1 " + 2V^V[) + 8^^) 

+L (r) J" dr'L c e 2Vo (W{M' Q + 4M (v(' + 2V£V{) + 8^^) (5.59) 
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The lower boundaries of the integrations are of course arbitrary, since their change can always 
be compensated by changing the integration constants m ,rhi. 

As r — > oo, the requirement that the perturbation V\ should maintain the asymptotics 
e 2V ~ 2cyr, requires that V\ must vanish at least as fast as 1/r. (The mode Vf, V™, and Vf 
all satisfy this requirement, while the dilation mode does not; thus the dilation mode must 
be absent, and we have Q = 0.) Given that M and L are linear in r as r — > oo, while Lq 
goes like 1/r, it follows that the leading asymptotics as r — > oo of the integrands is given by, 

L e 2Vo (2V[M' Q + 4M (v{' + 2V£V() + 8£ £i) ~ O(r ) 

L c e 2Vo (2V;M^ + 4M (v;' + 2V^V;)+8E E 1 ) ~ 0(r~ 2 ) (5.60) 

Hence, the second integral is convergent; the first integral behaves linearly as r — > oo so that 
the contributions of both integral terms behave as constants when r — > oo. Demanding that 
Cm be unchanged by the perturbation Et then requires that rhi = 0. 



5.7 Calculation of Ct at the critical point 



The coefficient Ct was already computed in section 4.3.3 by indirect matching methods. For 
completeness, we shall here also give its derivation directly from the perturbation theory 
solutions. The coefficient Ct is defined by the following relations as r — > 0, 



M(r) ~ — ar — 6tt 
N = C T e T 



(5.6i; 



while keeping ce constant. Actually, at the critical solution, we have a = 0. To determine 
the coefficient Ct, we need both the functions Mi and Ex, the latter because we need to 
make sure that when turning on iV at the horizon, we keep Ce unchanged. In fact, from 
formula (5.50), it is immediately obvious that to keep ce constant, we need to require 

C B 



ai = -no- 



CEo 

so that the field E\ = A\ relevant to the Et perturbation at constant ce is given by 



Ai(r) = A (r) 



-no / dr 

ceo 



A {r') V Lo L Q 



(5.62) 



(5.63) 



We begin by obtaining the r — > asymptotics of E\. Since we have Eo(r)/Ao(r) ~ k/r, the 
last two terms dominate and their integral produces the following leading asymptotics, 



A 1 (r) 



rioctk Ao{r) 



24V12(2A; - 1) 



O(r°)A (r) 



(5.64) 
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Now, the key ingredient we use here is that at the critical point, a = 0, so this singular 
leading asymptotics actually drops out. (Away from the critical point, we would have to 
compensate for this mode by an explicit subtraction of a translation mode.) As a result, the 
leading asymptotics of E\ at the critical point is just Ex(r) ~ r fc-1 . We are now in a position 
to evaluate the asymptotics of the integrands in (5.59) as r — > 0, taking V\ = Vf ~ Vq ~ r a , 
and Ei ~ Eq ~ r fc-1 . The terms in Mq and Mq scale as r a , while the term in scales as 



„2fc- 



; " 1 . As a result, the first integral in (5.59) is convergent as r — > 0, but in fact for k > 1/2, 
it vanishes faster than r as r — > 0, and thus never contributes to the linear behavior of 



M(r) which we need to extract following (5.61). Thus, the first integral is irrelevant to the 



calculation of Ct- The second integral in (5.59) is convergent as r — > as long as k > 1/2. 
Thus, we extract, 



1 



n jo 



dr' LqC 



2V 



(W{M' Q + 4M (V 1 " + 2V X) + 8£ £i 



(5.65) 



where V\ = Vf, and E\ = A\ is the field defined in (5.63). All dependence on n$ drops out of 
these formulas, and the remainder is expressed solely in terms of the functions of the T = 
critical solution times a factor of c% . 



6 Discussion 

In this paper we have restricted attention to the gravity side of the holographic duality, 
but it will of course be illuminating to further explore the field theory interpretation of the 
quantum critical point and its properties. In [11] we noted that the scaling properties of the 
entropy density agree with those near a quantum critical point described by the Hertz-Millis 
theory [3U |35l [36]. To develop this connection further it will be very useful to study finite 
frequency /wavelength perturbations around the critical point and to compute the associated 
dispersion relations. Based on our experience with finite temperature perturbations, it seems 
likely that this is analytically tractable. 

Even at the level of free field theory, there are possible origins for the emergence of a 
critical magnetic field. In the simplest interpretation, we assume that for B > B c the charge 
density in the field theory is due to fermions occupying the lowest Landau level. This is 
energetically favorable, as both higher fermionic Landau levels as well as bosonic levels have 
higher energy. As the magnetic field is decreased, or the charge density is increased, the 
maximal p% (the momentum parallel to the magnetic field vector) carried by states in the 
lowest Landau level increases, which increases their energy, E — \ps\. Eventually, the energy 
exceeds that of the higher fermionic levels and the lowest bosonic level. By dimensional 
analysis, this occurs when B 3 /p 2 is of order unity, with the precise number depending on the 
charge spectrum. It is possible that the critical point we see on the gravity side is in some 
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way connected to this, although free field theory clearly seems incapable of explaining the 
details, such as the s ~ T 1 / 3 scaling law. 

Another interesting question concerns the effective spacetime dimensionality of the fixed 
point theory, by which we mean the number of dimensions along which long range correlations 
are present. Several features seem to point to a 1 + 1 dimensional interpretation, especially 
the appearance of near horizon geometries that factorize as M 3 x R 2 . However, examining 
more closely our analysis of the low temperature thermodynamics, we see that this is not 
so clear cut. At nonzero charge density and magnetic field, and for vanishing P3 chemical 
potential corresponding to the vanishing of L at the horizon, there is no finite temperature 
solution taking the form of a three dimensional black hole times R 2 . Indeed, in Appendix |A| 
we write down the general factorized solution, and none obeys these conditions, except when 
k — 1. Our low temperature solutions always involve some excitation of V, corresponding 
to a fully five-dimensional solution. This seems to suggest that the degrees of freedom 
contributing to the low temperature entropy density live in more than 1 + 1 dimensions. To 
better understand this, it will again be very useful to have results for the spectrum of finite 
frequency /wavelength perturbations around the critical point. 
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A General solution with constant V 

In this appendix we derive the most general solution for which V is constant, and 6^0. By 
rescaling x\ and X2 we set V = 0. The reduced field equations are then 

Ml (NE + LP)' + 2kbP = 

M2 (LE + MP)' - 2kbE = 

El L" = AEP 

E2 M" = -AE 2 

E3 N" = -AP 2 

EA (L') 2 = M'N' + Ab 2 

CON MP 2 + 2LEP + NE 2 -6+^(L') 2 -^M'N' + b 2 = (A.l) 

We begin by eliminating the bilinears in E and P between equations El, E2, E3, 

(L") 2 - M"N" = (A.2) 
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If N' 7^ 0, we solve for M' in terms of V and N' using E4, 



M< = <^ (A.3) 



Differentiating once and substituting M" into (A.2), we obtain a factorized equation 



(N'L" - L'N" - 2bN") (N'L" - L'N" + 2bN") = (A.4) 

At least one of these factors must vanish. Without loss of generality, we choose this to be 
the first factor (the case of the second factor vanishing corresponding to b — >■ —b), so that 
N'L" — [L' + 2b) N" = 0. This means that the functions N' and V + 2b must be proportional 
to one another. Since we have assumed that N' ^ 0, we have in all generality, 

L' + 2b = aN' (A.5) 

for some constant a. As a result, we have L" = aN" which implies P(E + aP) = by 
equations El and E3. By an SL(2,R) transformation, any solution E + aP = is equivalent 
to a solution P = 0. UN' — 0, then we must also have P = 0. We conclude that all solutions 
with V constant and 6^0 are SL(2,R)-equivalent to a solution with P = 0. We shall now 
proceed to solve the general case where V is constant and P = 0. 

A.l The general solutions for P = 

Constant V and P = simplify El, E3 and CON to the following relations, 

L" = N" = NE 2 = Q-2b 2 (A.6) 

the first two of which are solved by 

L{r) = lo + hr 

N(r) = UQ + nxr (A. 7) 

for constant Iq, li, n , ri\. We now have the following three classes of solutions. 

1. When n\ = 0, and n = 0, we have b 2 = 3, with l\ = ±26, and the only non-trivial 
remaining equations are M2 and E2, which are solved by, 



E(r) = qr 



±k-l 



2q 2 

M(r) = m + m x r - , -r ±2k (A. 

k(2k =)= 1) 



for arbitrary constants q, m ,mi. This solution is equivalent to (2.19). 
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2. When ri\ = 0, and no 7^ 0, we rescale to no = 1, so that E = q is a constant obeying 
g 2 + 2b 2 = 6, with L = ±2b from E4, the constraint q(k =p 1) = from M2, and 



M(r) 



mo + mxr — 2q 2 r 2 



(A.9) 



for arbitrary constants mo and mi. This is equivalent to the k = ±1 solution of (2.23). 
3. When m ^ 0, we have M' constant by E4, E = by E2, b 2 = 3 by E4 and CON, and 

M(r) = m + m x r (A.10) 
with If — mini = 12. This solution is SL(2,R)-equivalent to AdS^, with a standard 



presentation given by li = y!2 and mi = ni = 0. 



B Solving for V perturbations on asymptotically AdS§ 

To solve for the perturbations V\, given that the perturbations N x and Pi are now known, 
we need to solve a coupled set of equations for V\ and fi, such as fV and E4'. Eliminating 
fi would result in a 4-th order differential equation in Vi, whose solvability is unclear. Here, 
instead, we shall take a less direct approach; it will result in a single third order differential 
equation for the full composite / e 2V ; this equation will then be linearized around the fo, Vq 
solution; two of its solutions will be known as the corresponding translation and dilation 
modes, and the third solution may then be constructed by quadrature only. 

Our starting point will be the following two independent equations for the full fields, 








fe 



2V 



24e 



2V 



(b.i; 



fiy'Y + f'V + J + 3e~ w - 6 + MP 2 + 2LEP + NE 2 



where f,g were defined in (2.11). The first equation coincides with (2.12), while the second 



is the original CON equation of (2.10). These two equations combined are equivalent to El, 
E4, and CON. 



B.I Derivation of a single third order differential equation 



Equation (2.16) allows us to express g in terms of /, /' and the combination A — ptu, which 



in turn may be evaluated with the help of the first integrals of (2.15). Carrying out these 



substitutions, we may re-express the second equation of (B.I) as follows, 



f{ y'Y + f'V + 



(/ 



l\2 



+ b 2 e~ AV - 6 + MP 1 + 2LEP + NE 2 



16/ 



4/ 



ll - /io^o + ^kb (2l AC + /i C 2 + vqJY 



-4V 



(B.2) 
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Next, we solve the first equation in (B.l) by parametrizing the most general solution with 
the help of a function X, such that, 



fe 



2V 



~2V 



2AX 
X" 



and eliminate / and V in favor of X in (B.2). This gives the following equation, 

- 3(XX"') 2 + 6XX'X"X"' + (X'X" f + 2XX" - 4X(X") 3 
where R is defined by 



R 



(B-3) 



(B-4) 



R 



— — (MP 2 + 2LEP + NE 2 )fe 8V 
36 



144 



AV 



(B.5) 



In the perturbative expansion of this equation, R will act as source terms. For this reason, 
we have refrained from recasting its / and V dependence in terms of X. 

The unperturbed solution in terms of L , M , Vb, A has N = Co = Po = 0. Its boundary 
conditions A (0) = L (0) = M(0) = force the values of the first integrals \o = fiQ = u = 0, 
so that R = 0, and X obeys 



- 3(X X-) 2 + GXoX'XX' + {X'^f + 2X Q X>> - 4X (^') 3 = 



(B.6) 



First order perturbation theory around this solution was organized in (5.35), and we set 



X 
R 



X + eX 1 
R + en Ri 



(B.7) 



By inspection of (2.15), the constants Ao,yUo and u are all at least of order e (since they 
vanished for the unperturbed solution). Thus, Ri may be readily evaluated and we find, 

(B.8) 



Rx = --^(2L £ Pi + iVi^)/oe 8 ^ + ^Ale^ 
6bn 6bn 



The constant z/ is readily deduced using the last equation of (2.15), the expression for Ni 

(B.9) 



in (5.39), as well as (5.33), and we find, 



-no 



Finally, we linearize also the left side of (B.4), and organize the result as follows, 



a 3 X'" + a 2 X" + a x X[ + a X 1 = hoRt 



(B.10) 
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where the coefficients are found to be, 

a 3 = — 6XqXq" + 6X X' Xq 

a 2 = 6X X^" + 2(X^) 2 X;' + 2X -12X (Xj) 2 

01 = qx x , ( ;x^ + 2x' (x' ( ;) 2 

a = -6XoTO 2 + 6XKXj' + 2X;'-4(X;') 3 (B.ll) 
It is this equation for X\ than remains to be solved. 



B.2 Constructing the general solution for X\ 



Equation (B.6) for the unperturbed solution is invariant under translations and dilations in r. 
As a result, the homogeneous part of (B.10) manifestly admits two solutions: the collective 
coordinate modes X\ and Xf associated respectively with translations and dilations of the 
unperturbed solution X . We know their form explicitly in terms of X , 

X\ = X' Xj = tXq — 2Xq 



(B.12) 



The fact that X\ solves the homogeneous part of (B.10) allows one to reduce the order of 
the differential equation by one. To do so, we introduce a function Y which satisfies, 



X' Y 



b,Y"' 



b 2 Y" + b x Y' = h Q R x 



(B.13) 



where 



h 
61 



o 3 Xq 

CL 2 Xq + 3a3XQ 

ai X' + 2a 2 Xl + 3a 3 X;" 



(B.14) 



Since the homogeneous part of equation (B.10) admits also a second solution, namely Xf, it 
follows that the homogeneous part of the equation in Y must admit the solution, 



Yn 



rX' - 2X 



Introducing now a function Z which satisfies, 

Y' = Y^Z Z'%Y^ + Z\2bX + b 2 Y^) 

The homogeneous part of this equation reduces to, 

Z'n _ Y n _ Xn a 2 



n Ri 



° 7 + 2^ + 3 



'X'o 



+ 







«3 



(B.15) 



(B.16) 



(B.17) 
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The quantity 02/03 may be evaluated in terms of the functions Lq and Vo, and may be 
reduced to the following simple expression, 

- = -<-^-§ (B.18) 

Putting all this together allows us to compute, 

L L' e AVo 

z '° = (R19) 

Finally, it remains to find a particular solution to the inhomogeneous equation for Z of 
(B.16). To do so, we introduce a function Q which satisfies, 

Z' = n»Z' Q x Q\ = * l (B.20) 

a 3 A r Z 

Clearly, this equation can be solved by quadratures only 

We are now ready to present the general solution for X 1; in terms of its three linearly 
independent homogeneous solutions Xf, Xf, and the third solution which we shall denote 
X™, as well as a particular solution Xf sourced by R\. They are given explicitly as follows, 

X\(r) = X' (r) 
X*(r) = X f (r)Y Q (r) 

X-{t) = X' (r) ^ dr%(r')Z (r') 

X{{r) = X' (r) fdr%{r') [' dr'%(r")Qi(0 (B.21) 

JO Joo 

where we have defined Qi to vanish at r = 0, 

Qi(r) = f dr'Q'^r') (B.22) 
Jo 

The general solution for X% is then given by, 

X t = CtXl + CdXt + CnX? + n Xf (B.23) 

The scripts t,d,n,p refer respectively to translations, dilations, non-local, and particular 
solution, and (t, (d, and ( n are integration constants. The normalizations are such that 
Qi(0) = 0, K (0) = 0, and Z (oo) = 0. 
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B.3 Asymptotic behavior of X\ and V\ 

The effects on the fields of the translation and dilation modes are well-known, and have 



already been spelled out in (5.38). To derive the asymptotics of the Xf and Xf modes, 
and their corresponding effects on the other fields, we express the modes in terms of the 
unperturbed solution and the fields N\ and P\. One finds, 



1 



X> = -L (L> + L V')e 2V ° 



y 



Z' 



12 - (L' + L V>) 2 

(L' + L vzy 

(12) 3 L' L' + LqK 
L 2 e 2V o [12 - (L' + L V{) 



[B.24) 



result, the equation for Q 1 becomes, 

12 - (L' + L Vq) 2 



Q'i 
Ri 



3Lq(L' 



2 p 6V 



Ri 



1 



36fi 



^(2L E P 1 + N l E 2 )L 2 e 



2„8V 



— A 2 e m 
36 



(B.25) 

The r — > oo asymptotics of X 1 may be evaluated using the asymptotics of V , L , and of 



Ni,Pi from (5.40) and (5.45) respectively, and we find 



X{ 

xf 



t °V r 2 

2 



6 '' 



24r 



where Ji(fc) is defined by the integral 

Ji{k) = 



dr'Y^r') 



MQXl + QiWX? 



dr" Z^{r")Q x (j") 



(B.26) 



(B.27) 



which is convergent for all k > 1/2. Clearly, the boundary conditions as r — > oo require the 
dilation mode to be absent, so that Q = 0, but allows for the three remaining modes. 

The r — » asymptotics of X\ may be evaluated similarly. Actually, we shall be interested 
only in contributions to X± that produce functional dependences in V\ of the form r u and 



-l-cr ; 



in this limit, since these are the only modes that enter in the overlap region of (5.16) 



Since we have 2Vie 2l/o = X", the corresponding modes in X\ behave as r a+2 and r 1 a . Using 



the r — > asymptotics of Lq,Vq, and Ni,P\ from (5.40) and (5.45) respectively, we find 



X\ ~r 

Xf ~ ^Zir^ 
1 a + 1 



Xf 
Xf 



6 r 1 



(j(2(7 + l)(3(7-2) 

Mk)x* 



[B.28) 
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where J^k) is defined by the integral 

POD 

J 2 (k) = - / drZ^Q^r) (B.29) 
Jo 
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